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High Performance Time-Frequency 
Distributions for Practical Applications 

Boualem Boashash 
Victor Sucic 

ABSTRACT This chapter presents in three interrelated sections the key concepts 
and techniques needed to design and use high performance time-frequency distribu- 
tions (TFDs) in real-world practical applications. 

Section 6.1 first presents, in a heuristic approach, the core concepts forming the 
field of time-frequency signal processing, incorporating recent developments, such 
as the design of high resolution quadratic TFDs for multicomponent signal analysis. 

Section 6.2 outlines methods of assessment of the performance of time-frequency 
techniques, in terms of the resolution performance of TFDs in separating closely 
spaced components in the time-frequency domain. A performance measure is defined 
using key attributes of TFDs, such as the components’ mainlobes and sidelobes, and 
cross-terms. This method of assessment of TFDs performance has led to improve- 
ments in designing high resolution quadratic TFD for time-frequency analysis of 
multicomponent signals. 

Section 6.3 presents a methodology for selecting the optimal TFD for a given 
real-life signal under application-specific constraints. The methodology, based on 
the performance measure, allows for emphasis of signal features in specific regions 
of interest in the time-frequency domain. 



6.1 The core concepts of time-frequency signal 
processing 

6.1.1 Rationale for using time-frequency distributions 

In order to provide more insight into the nature of nonstationary signals, a new field 
of science and engineering has emerged: the field of joint time-frequency signal pro- 
cessing (TFSP). 

The introduction of TFSP has led to new tools to represent and characterize the 
time-varying contents of nonstationary signals using time-frequency distributions 
(TFDs), the most popular ones belonging to the class of quadratic distributions (see 
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Eq. (6.1 .20)). By distributing the signal energy over the time-frequency plane, TFDs I 
provide the analyst with information unavailable from the signal time-domain repre- 
sentation or its frequency-domain representation. This includes the number of com- 
ponents present in the signal, the time durations and frequency bands over which 
these components are defined, the components’ relative amplitudes, phase informa- j 
tion, and the instantaneous frequency (IF) laws that components follow in the time- 
frequency plane. 

The essential characteristic of TFSP is that it comprises a set of signal process- 
ing methods, techniques, and algorithms in which the two natural variables time, t, 
and frequency, /, are used concurrently. This contrasts with the traditional signal 
processing methods in which time and frequency variables are used exclusively and 
independently. Nature shows us in our daily experiences that the two variables, f and 
/, are usually simultaneously present in signals (e.g., the speech of a person, the 
song of a bird, or music played on the radio). Such signals are called “nonstationary" 
signals because their spectral characteristics vary with time. This chapter considers 
signals that have finite duration T and nearly finite bandwidth B; these signals are 
often referred to as asymptotic, with the degree of asymptoticity expressed by the 
BT product [3]. TFSP is designed to deal effectively with such signals by allowing 
a more detailed analysis and processing. 

In addition to existing methods of TFSP, special-purpose TFDs could be defined 
to suit the particular class of signals under investigation. This chapter presents a 
simple approach for TFD design. For this purpose, we first revisit the core concepts 
of TFSP from a practical point of view (Section 6.1), and investigate the performance 
criteria for TFDs from a user’s point of view. The performance of TFDs [1,2] is 
assessed in terms of concentration and resolution , and an optimization procedure is 
used to select the optimal parameter values of TFDs (Section 6.2). Finally, we define 
a methodology for selecting the most suitable TFD in a given practical situation 
(Section 6.3). 

6.1.2 Limitations of traditional signal representations 

The spectrum of a signal (deterministic or random) gives no indication as to how the 
frequency content of a signal changes with time — information that is important when 
one deals with frequency modulated (FM) signals or other kinds of nonstationary 
signals. This frequency variation often contains crucial information about the signal 
and process studied in applications. 

The limitation of ‘classical” spectral representations is better illustrated by the 
fact that we can find totally different signals (related to different physical phenom- 
ena), s a (t) and Sb(t), which yet have the same “spectra” (that is, magnitude spectra). 

The following example shows the inherent limitation of conventional spectral analy- 
sis. 

Example 6.1. Consider the two signals, s a (t) (finite duration linear FM) and s b (t) 
(infinite duration sine function), defined as 
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(a) Time domain 



(b) Magnitude spectrum 



(c) Phase spectrum 



Figure 6.1. Time-domain and frequency-domain representations of the finite duration linear 
FM signal s a (t) defined by Eq. (6.1.1). 




0<t<T 

elsewhere 



( 6 . 1 . 1 ) 

( 6 . 1 . 2 ) 



where /q is the start frequency, a = B/T represents the rate of the frequency change, 
and f c ~ fo + B/2, 

The signal s a (t) is nonstationary: it is a linear FM signal commonly used in 
radar and seismic applications (it is analogous to a musical note with a steadily ris- 
ing pitch). Equation (6.1.1) indicates that the signal s a (t) is a cosine function in 
the interval (0, T), and zero outside this interval. Figures 6.1 and 6.2 show the time 
representation and the frequency representation of the signals s a (t) and s&(i), re- 
spectively. 
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(b) Magnitude spectrum 




(c) Phase spectrum 



FIGURE 6.2. Time-domain and frequency-domain representations of the infinite duration sine 
signal s;,(t) defined by Eq. (6.1.2). 
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Although the signals s a (t) and Sb(t) are fundamentally different, they have the 
same magnitude squared spectrum (Figures 6.1(b) and 6.2(b)). The information that 
allows one to discriminate between them is contained in their phase spectra (Fig- 
ures 6.1(c) and 6.2(c)), which is lost when we form the square modulus of the Fourier 
transform (FT), The frequency in the signal s a (t) is steadily rising with time, a fact 
not revealed by the signal spectrum displayed in Figure 6.1(b), which only shows a 
constant amplitude broadband spectrum. 

It follows then that the magnitude spectrum by itself is insufficient (and hence 
inadequate) for representing signals in a way useful for precise characterization and 
identification. This serves as a part of the motivation for devising a more “sophisti- 
cated” and “practical” nonstationary signal analysis tool, which preserves “all” the 
information of the signal and discriminates signals in a better way, using a single 
complete representation instead of attempting to interpret signals’ magnitude and 
phase spectra separately. 

6.1.3 Positive frequencies and the analytic signal 

An important property of the FT is that for s(t) real, its FT, S(f), is complex and 
has Hermitian symmetry (S(f) = S*(-f )); i.e., its real part and magnitude are even 
symmetric, and its imaginary part and phase are odd symmetric. The information 
contained in the negative frequencies is therefore a duplication of the information 
contained in the positive frequencies. As a frequency represents the number of oscil- 
lations of a signal observed in a second, a frequency is expected to be positive. The 
negative frequencies appear in the spectrum of a real signal as a consequence of the 
mathematical model of the FT. In practice, a signal analyst would prefer to deal only 
with positive frequencies. This is achieved by introducing a complex signal called 
the analytic signal. 

Definition 6.1. An analytic signal z(t) is a complex signal that contains only positive 
frequencies, such that its FT, Z(f), satisfies 

z (f) = 0, for / < 0. (6.1,3) 

A real signal x(t) cannot have only positive frequencies, as the FT of x(t) has a 
Hermitian symmetry, which contradicts the hypothesis of X(f) - 0 for / < 0. 



Theorem 6.1. A complex signal z{t) = x{t)+j y{t) is analytic if and only if its real 
and imaginary parts are related as follows: 

y(t)= x {t)*~=n{x(t)}, 

where * denotes the convolution in time t and H{f is the Hilbert transform: 



(6.1.4) 



6. High Performance TFDs for Practical Applications 139 
= p.v. | J°° dr j (6.1.5) 

with p,v- being the principal value of Cauchy’s integral. 



Theorem 6.2. // s(f) « a rcrd signal expressed as s(t) = a(t ) cos then un- 
der the conditions outlined in [4], and using the fact that H{a(t ) cos <?5(t)} = 
a(t)H{cos 4>{t)}, we can construct an analytic signal z(t ) which corresponds to 
s(t) by replacing cos <f>(t) by e j4 ‘ < ~ t '>; that is, z(t) = a(t) and s(t ) is the real 
part of 2 (f). 

This theorem is important for practical applications because it shows a relation 
of equivalence between the real signal s(t) and its analytic associate z(t). 

6.1.4 Joint time-frequency representations 

6. 1 .4. 1 Uncovering hidden information using time-frequency representations 

Revealing the time and frequency dependence of a signal, such as the signal s a (t) 
in Example 6.1, is achieved by using a joint time-frequency representation. In Fig- 
ure 6.3, the Wigner-Ville distribution (Eq. (6.1.16)) is used to represent s a (t), as 
it provides the optimal joint time-frequency representation for a chirp (linear FM) 
signal [3], 

In Figure 6.3, the signal start and stop times are clearly identifiable from its time- 
frequency representation, as is the time variation of the frequency content of the 
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Figure 6.3. Time-frequency representation (Wigner-Ville distribution) of the linear FM sig- 
nal s G ({) defined by Eq. (6.1.1). The signal’s time domain representation is given on the left, 
and its magnitude spectrum on the bottom. 
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(a) Modified B distribution of the first (b) Modified B distribution of the 
signal second signal 

Figure 6.4. Time-frequency representation (Modified B distribution) of two different three- 
component signals. 



signal. This information cannot be retrieved solely from either the signal instanta- 
neous power |s(i)| 2 or its spectrum |S(/)| 2 representation. The ease of interpreting 
plots such as the plot in Figure 6.3 makes the concept of a joint time-frequency signal 
representation attractive. 

Example 6.2. Figure 6.4 illustrates another case of two different signals having sim- 
ilar magnitude spectra. Both signals contain three linear FM components. The differ- 
ences in the time intervals, frequency bands, and the FM laws that characterise the 
signals components are not shown clearly in the t domain and the / domain, but do 
appear clearly in the signals’ (t , /) representation (here the Modified B distribution 
(MBD) [5]). 

Example 6.2 shows that another advantage of using a joint time-frequency repre- 
sentation of signals is that it reveals whether the signal is monocomponent or multi- 
component [3,6], a fact that cannot be easily obtained from the signal time-domain 
or frequency-domain representation. 



6. 1 .4.2 What is a time-frequency representation? 



TFSP is a natural extension of both time-domain and frequency-domain process- 
ing, that involves representing signals in a complete space that can display “all” the 
signal information in a more accessible way [6], Such a representation is intended 
° provide & distribution of signal energy versus both time and frequency simultane- 
ously. For this reason, the representation is commonly called a TFD, and is denoted 



auH fi C °Ti Pt “ y rdated t0 joint time ~ frec l ue nc y representation is that of IF 
and time delay. The IF corresponds to the frequency of a sine wave that locally (at a 
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given time) fits the signal under analysis. The time delay is a measure of the “order 
of arrival” of the frequencies. 

The TFD is expected to visually exhibit in the (t, /) domain the time-frequency 
law of each signal component, thereby making the estimation of their IFs and time 
delays easier, and should also provide additional information about components’ rel- 
ative amplitudes, and the spectral spread of the components around their IFs (the 
spread is known as the instantaneous bandwidth [3, 6, 7]). For example, the ridges in 
Figure 6.4 indicate that the three signal components have equal amplitudes, and the 
peaks of these ridges reveal the linear FM laws of the components. All components 
have the same instantaneous bandwidth around their IFs, as indicated by the same 
concentration the ridges attain about their peaks. 

6.1. 4.2.1 Instantaneous frequency and time delay 

The instantaneous frequency of a monocomponent signal is a measure of the local- 
ization in time of the individual frequency components of the signal [4], 

The IF, of a monocomponent analytic signal z(t) — a(t)e is given by 

(61 ' 6) 

The IF of a monocomponent real signal s(t) = a(t) cos <f>(t) is defined as the IF 
of the analytic signal z(t) corresponding to s(t). 



The dual of the IF concept in the frequency domain is called the “time delay.” 
The time delay, Td{f), of a monocomponent analytic signal z(t) is defined as 



r d {f) = - 



1 dd(f) 
df ' 



(6.1.7) 



where 



= Z(f) = ( 6 - 1 - 8 ) 



is the FT of z(t). 

The time delay of a monocomponent real signal s(t) is defined as the time delay 
of the analytic signal z(t) corresponding to s(t). 

The order of appearance of each time-varying frequency component is the time 
delay. The global order of appearance of the frequencies is called the group delay 
(a mean value of individual time delays). The group delay equivalent for IFs is the 
mean IF. 

The IF and the time delay describe the “internal organization” of the signal. The 
use of time-domain representations or frequency -domain representations leads one 
to effectively neglect this information, resulting in a scrambling of the information 
contained in the signal, as illustrated in Examples 6.1 and 6.2. 
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6.1.4.2.2 Physical interpretation ofTFDs 



Most TFDs used for a practical time-frequency signal analysis are not necessarily 
positive, so they do not represent an instantaneous energy spectral density at time f 
and frequency /. For example, the Page distribution is defined as the time derivative 
of the “running spectrum” (spectrum of the signal from — oo to time t), and hence 
can take both positive and negative values [8]: 



*<•■/>- s 



z(u)e~ j2 * fu du 



To relate p z ( t , /) to the physical quantities used in practical experimentation, we 
interpret p z (t,f) as a measure of energy flow through the spectral window 
(/ - A//2,/ + A//2) during the time interval (t - At/2, t + At/2). The sig- 
nal energy localized in the time-frequency region (At, A/) is then given by [3] 



pt+At/2 pf+Af/2 

EAt,Af- / p z (t,f)dfdt. (6.1,9) 

Jt-At/2 Jf-Af/2 

The spectral window should be chosen large enough so that the product A tAf sat- 
isfies the Heisenberg uncertainty relation [9]: 



AfA/> l/(47r). (6.1,10) 

With this interpretation, negative values of p z (t, f) are then accounted for, and posi- 
tivity can be removed as a requirement. This is especially important since positivity 
is also incompatible with the requirement that the signal IF is the first moment of its 
TFD with respect to frequency [3], 

6.1.4.2.3 Desirable properties and criteria for a TFD 



The TFD is expected to satisfy a certain number of properties that are intuitively 
desirable for a practical analysis. It was reported in [3] that a TFD should satisfy 
the marginals; i.e., a TFD should reduce to the spectrum and instantaneous power by 
integrating over t, respectively / , Furthermore, a TFD is expected to have the IF as 
its first moment with respect to frequency. These strict constraints on the TFD design 
led to the terminology of Cohen’s class [10], 

However, our approach is different and more in line with actual usage and prac- 
tice. We first note that many popular TFDs (e.g., the spectrogram) do not satisfy the 
marginals and the IF moment condition. Yet the spectrogram has been a valuable tool 
m many practical applications, suggesting that the time and the frequency marginal 
and the IF moment constraints may not be strictly needed in practice. What is more 
important in most practical applications is to maximize the energy concentration 

about the IF for monocomponent signals and improve the resolution for multicom- 
ponent signals. 

We have found from practical experience that a time-frequency representation of 
signal z(t), p z ( t , /), needs to satisfy the following properties, so that it can be useful 



forpra 

• PI (J 
time-f 



» P2 ( 
plane, 
of the 



Ezn i 
Note i 
(Eq.( 

• P3 i 
comp 



Forn 

ponei 

• P4 
expei 
comp 
signs 
pone 

. P5 
analy 
plani 
nate 
the 1 
It is 
com] 
estin 

1 

Pl-1 
and 
TFE 
are i 



6. High Performance TFDs for Practical Applications 143 



for practical purposes and not just for theoretical interest [2, Part I]. 



• PI ( Reality and energy ): The TFD should be real, and its integration over the whole 
time-frequency domain should be equal to the energy of the signal z(t), E z : 




Pz{tJ)dtdf = E z . 



( 6 . 1 , 11 ) 



* P2 ( Distribution of energy): The signal energy in a certain region R in the (t, f) 
plane, E ZR , should be obtained by integrating p z (t, f) over the boundaries (At , A/) 
of the region R: 




dt = E zr . 



(6.1.12) 



E Zfl is a portion of signal energy in the frequency band A / and time interval At. 
Note that A / and At need to be selected in such a way that the uncertainty principle 
(Eq. (6.1.10)) is satisfied. Property P2 is equivalent to Eq. (6.1.9). 



• P3 (IF peak property): The peak of the time-frequency representation of a mono- 
component FM signal with respect to frequency should reflect the IF of the signal: 



dpzjtj) 

df 



= 0 . 

/=/«(t) 



(6.1.13) 



For multicomponent signals, the same property should apply to the individual com- 
ponents. 



• P4 ( Concentration and resolution): The TFD of a monocomponent FM signal is 
expected to have a good energy concentration around the signal IF law. For a multi- 
component FM signal, a TFD is expected to provide a good (t , f ) resolution of the 
signal components. This requires a good energy concentration for each of the com- 
ponents and the suppression of any undesirable artifacts. 

♦ P5 ( Robustness to noise): When a signal embedded in additive white noise is 
analyzed in the joint time-frequency plane, the noise becomes evenly spread in the 
plane [3], For moderate to high signal-to-noise ratios, the signal components domi- 
nate over noise in the time-frequency domain, so their IF laws can be estimated from 
the TFD’s dominant peaks (Property P3), provided the TFD satisfies Property P4. 
It is desirable that such obtained IF estimates are as close as possible to the signal 
components’ true IF laws, i.e., that they are minimum bias and minimum variance IF 
estimates. 



Naturally, the following questions arise: Is there a TFD that meets specifications 
P1-P5? If yes, how do we construct it? What are the significant signal characteristics 
and parameters that impact on the design of such a TFD? How do these relate to the 
TFD itself? How do we obtain them from the actual designed TFD? These questions 
are answered in the following sections of this chapter. 
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6.1.5 Heuristic formulation of the class of quadratic TFDs 

6. 1 .5. 1 Time-varying spectrum and the Wigner-Ville distribution 



Considering the basic definition of the power spectrum density (PSD), let us deter- 
mine why time information seems to disappear when we take the PSD, and how it 
can be restored. 

Let us consider a complex random process Z(t) with realizations z(t, e), where e 
represents the ensemble index identifying each realization. Using context, rather than 
explicit rigorous mathematical formulation, to improve clarity we simply replace 
z(t, e) by z(t). Thus, e is implicit when we say that z(t) is random. 

The Wiener-Khintchine theorem states that for a stationary signal, the signal 
PSD equals the FT of its autocorrelation function [1 1]. By extension to nonstationary 
random signals, the time- varying PSD, S z (t, /), is then defined as the FT of the time- 
varying autocorrelation function, R z (t , r). 

The time-varying autocorrelation function of z(t ) is defined in symmetric form 
as [3] 



R z {t,r) = £{z{t + r/2)z*(t -r/2)} 

= £{K z (t )T )}, (6.1.14) 

where K z (t, r) = z(t + r/2)z*{t - r/2) is the signal kernel and £{•} defines the 
expected value operator. 

Therefore, the signal time-varying PSD S z (t, f) is given as 



&(*,/)= {R z {t,r)} 

T~>f 

= T {£{K z (t,r)}} 

T—*S 

= £ { T Zf 

= £{W z (tJ)}, (6.1.15) 

where the interchange of £{■} and FT is made under assumptions verified by most 
practical signals and asymptotic signals [3]; and W z (t,f) denotes the Wigner-Ville 
distribution (WVD) of z(t), expressed as 



W z (t,f)= R{K z (t,r)} 

T-*f 



= R {z(t + r/2)z*(t — r/2)} 
r— ♦/ 



/ oo 

z{t + r/2)z*[t - r/2) e~ j2nfT dr. (6.1.16) 

-OO 



The problem of estimating the time- varying PSD of a random process z(t) thus re- 
duces to averaging the WVD of the process over e. If only one realization of the pro- 
cess is available, assuming the process is locally ergodic [1 1] over a time window, an 
estimate of the process time-varying PSD is obtained by smoothing the WVD over 
the window of local ergodicity [3], 
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6 . 1.5.2 Time-varying spectrum estimates and quadratic TFDs 

If z{ t) is deterministic, from Eq. (6.1.15), S z (t, f ) reduces to W z (t, /): 

s,( t, f ) = £{W z (t , /)} = W z (t, /). (6.1.17) 

The signals we consider in this chapter are asymptotic (most real-life signals), having 
a finite duration T, and a finite bandwidth B [3], The finite duration of a signal is 
obtained by windowing the signal with a finite duration time function g(t), hence 
convolving S z (t , /) in tire frequency domain with G(f) = T 

t >/ 

The finite bandwidth restriction is met by windowing the signal in the frequency 
domain with a band-limited function H(f), hence convolving S z (t, /) in the time 
domain with h(t) — T 

By combining the separate windowing effects of g(t) in time and H(f) in fre- 
quency, the estimate of the signal tune- varying PSD can be defined as 

S z (t, f) = G(f) */ W z (t, f) * t h(t ), (6.1.18) 

where G(f) and hit) are even functions (such as those traditionally used in spectral 
analysis and digital filter design), and and */, respectively, denote the convolution 
in time and convolution in frequency. 

The preceding formulation for S z ( t , f) was introduced step by step to illustrate 
the two-dimensional (2D) convolution that is inherent to most real-life signals. This 
formulation, however, corresponds to a special case where the 2D windowing in t 
and / is separable. 

In the general case, we need to introduce an even function 7 (t, /), which may or 
may not be separable, that reflects the signal overall duration-bandwidth limitations 
in both time and frequency. This leads to the following general formulation of time- 
frequency representations p z (t, f) of signal z(t) [2, Part I]: 

S z (t, /) = p z {tj) = W z (t,f)* t * f 7 (*,/), (6-1.19) 

where 7 (t, f) is an even function that defines the TFD p z (t, f) and its properties, 
and denotes the convolution in both time and frequency. 

By expanding W z ( t , /) and the double convolution in Eq. (6. 1. 19), we obtain the 
following general quadratic form of TFDs [2, Part I]: 



Pz{t, /) 



/ CO poo poo 
-00 J — 00 </ — 



j2nu(u 



— OO J — OO </ — OO 



j2irfT dududr , 

( 6 . 1 . 20 ) 



where cj(t, i v), known as the TFD kernel filter, is the 2D FT of 7 (t, /). 

The function cj(t, v) is analogous to the windows used in classical spectral anal- 
ysis, and it determines how the signal energy is distributed in time and frequency. 
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By appropriately choosing g(r t v) in Eq. (6.1 .20), we can obtain most of the popular 
time-frequency representations of z(t), such as those defined in [2, Part I], 

Eq. (6.1.20) differs from Cohen’s class [7, 10] by the sign in the first exponential. 
This allows for fewer restrictions on the choice of kernel filter, as discussed in Sub- 
section 6.2.3, and the interpretation of p z (t, f) as a filtered version of the WVD, as 
explained in Subsection 6.I.5.3. 

6. 1 .5.3 Time, lag, frequency, Doppler domains, and the ambiguity function 

We have defined in Eq. (6.1.1 6) the WVD of signal z{t) as the FT of the signal kernel 
K 2 (t,T ): 

W z (t,f) = (6.1.21) 

The representation W z (t, f) is a function of two variables, t and /. 

The time-averaged autocorrelation function of signal z(t ) may be defined as 

/ OO 

oo z ^ t+ ^ z *( t ~^) dt - ( 6 . 1 . 22 ) 

i? z (r) describes the similarity between the signal and its time-delayed copies. It is 
obtained by taking the integral over time of the signal kernel K z (£, t). 
r A not h er quantity that relates to K z (t,r) is the symmetric ambiguity function 



/ OO 

K z (t, t) e-V^dt 

•OO 

/ OO 

^ Z(t + T -)z*{t - T -) e-^dt. (6.1.23) 

U3) f n£S the basic radar e ^ uation obtained by correlating a signal 

Thus I? “ timC by kg T and Shifted in frequency by Doppler (a 

concern ’ W f “ I £lement “ the formuIa tion of many important TFSP 

rZcv ?™fuT aCy t J) ’ time ‘ lag (t ’ r) ’ la 8' D °PPfrr (r,%, and Doppler- 

means of he’pT and 10ns ° btamed fr0m the si 8 nal kernel K x (£, r) by 
means of the FT and the inverse FT [2, Part I], This is illustrated in Figure 6 5 

of s .gn um iriha 6 ' 5 ' the 2DFT<indicatedb y ,W0 vertical allows) of the WVD 
gnal equals the symmetric ambiguity function A z (r, v) of z(t) [2, Part I]: 

W ^ t >f) t f^r A z (t, v). (6.1.24) 

the in J FSP bECaUSe * indiCates the dual of 

domain). The It f ) domain ren £ ^ ° Ppler domain ( also called the ambiguity 
actual frequency, ancf the^T- *** “ 3 fUDCtion ° f actual ‘-e and 

shifts and frequency shifts. ’ represents the signal as a function of time 
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W z (tJ) 




Az{t,v) 



Figure 6.5. Different dual domain representations of signal z(t), obtained from the signal 
kernel K x ( t , r). 



Given that a double convolution in (t, /) results, by 2D FT, in a multiplication 
in the ambiguity domain, A z (t, u) becomes a windowed version of A z {t, v). This 
leads to the interpretation of TFD design, defined by Eq. (6.1.19), as a 2D filtering 
procedure in the ambiguity domain. 

Equation (6.1.20) can therefore be rewritten as the 2D FT of the symmetric am- 
biguity function A z (t, v) filtered by the TFD kernel filter g(r, v) [2]: 

Az(t, u) g(r, v) — A z (r, v) p z (t,f)=f f g(r,u)A z (r,i/)e :,2lT< -'' t ~ fT) dv dr. 

J — OO J —GO 

(6.1.25) 

Choosing the kernel filter g(r, v) in Eq. (6.1.25) most relevant to an application 
results in a specific TFD. For an all-pass filter, g(r, v) — 1, p z (t, f ) reduces to the 
WVD. For g(r, v) chosen to be the ambiguity function of a time analysis window 
w(t), p z (t, f) corresponds to the spectrogram with window w(t) [2,3]. 

6. 1.5.4 Quadratic TFDs, multicomponent signals, and cross-terms reduction 

Equation (6.1.20) defines TFDs that are quadratic (or bilinear) in the signal z(t). 
This implies that if z(t) includes two components z\(t) and 22 (f), then its quadratic 
formulation will not only include these two components but also additional compo- 
nents corresponding to their cross product zi (t) 22 (f)- These additional components 
are often called cross-terms and are considered as “artifacts” or “ghosts” appearing 
unexpectedly in the (t, /) representation. In Figure 6.6 (the WVD of two linear FMs), 
the “ghost” component is located halfway between the two signal components, and 
has an (oscillating) amplitude that exhibits large positive and negative values [3]. 
An effect similar to cross-terms appearance in (f, /) representations occurs when we 








148 B. Boashash and V. Sucic 




Figure 6.6. WVD of a signal composed of two linear FMs. The middle “component,” ex- 
hibiting large positive and negative amplitudes, is known as the cross-term. 



take the spectrum of zi(t) + Z 2 {t)\ we obtain cross-spectral components that are zero 
only when z\(t) and z 2 (t) do not overlap in frequency. 

Thus, the introduction of either noise or some other deterministic components 
in z(t) introduces cross-terms in its representations. In some applications the cross- 
terms may be useful, as they provide additional features that can be used for signal 
identification and recognition [6, Chapter 18], However, in most cases, they are con- 
sidered as undesirable interference terms that distort the reading of the representa- 
tion. So it is generally desired to design TFDs that suppress them “best.” 

6.1. 5,4.1 Location of cross-terms 

It was shown that a signal mapped by A z (r, v) into the lag-Doppler domain always 
traverses the origin of that plane, while the cross-terms, having oscillating amplitude 
in the time-frequency domain, are located away from the origin in the lag-Doppler 
plane, the distance being directly proportional to the time and frequency separation of 
the signal components [2, Articles 4.2 and 5.1], [12], This is illustrated in Figure 6.7 
for a two-component signal whose (Gaussian-like) components are centered in the 
time-frequency domain at (h , f x ) and (t 2 , f 2 ), respectively. 

Since the WVD is related to the ambiguity function by a 2D FT (Eq. (6.1 .25)), the 
simplest way to reduce the cross-terms of the WVD would be by filtering them out 
in the ambiguity domain, before taking a two-dimensional FT to return to the (it, /) 
domain. The two-dimensional filter g(r, v) needs to be chosen such that it “passes” 
the legion of the ambiguity plane close to the origin (where the signal components 
are located), and at the same time attenuates the rest of the plane. Note that any 
resulting truncation of the signal tenns by the kernel filter g(r, v) would result in the 
signal components spreading in the (t,/) domain, and so in the loss of their time- 
frequency resolution. So, a compromise is needed when defining the kernel filter in 
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Figure 6.7. Location of signal terms and cross-terms in the (t, /), left, and (r, v), right, 
domain. 



the ambiguity domain (and setting its parameters) such that the best possible trade- 
off is achieved between time-frequency resolution and cross-terms suppression. 

6. 1.5. 4. 2 Reduced interference TFDs 

, Many TFDs have been designed for the purpose of cross-terms suppression, the most 
} popular among them belonging to the class of reduced interference distributions 
(RIDs) [6, 13]. One of the first RIDs was the Choi-Williams distribution (CWD), 
whose g (r, u) is a 2D Gaussian function centered around the origin in the ambiguity 
plane, with its spread controlled by a parameter a [13]. 

Another, more recent, reduced interference TFD is the MBD [5], whose time-lag 
kernel filter is defined as 

= (6 ' L26) 

where k — F (2/5) / (2 2t3 ~ 1 T 2 (p)) is the normalizing factor (F(-) stands for the 
gamma function), and /3 is a real, positive number that controls the sharpness of the 
cutoff of the kernel filter in the ambiguity domain, and so also controls the trade-off 
between time-frequency resolution and cross-terms suppression. 

The MBD was found to be the closest to the ideal compromise; it is almost cross- 
terms free and has high components’ resolution in the time-frequency plane [5], In 
addition, the MBD allows an efficient IF estimation for multicomponent signals [5]. 
The original B distribution (BD) [14] also performs well, but it does not meet the 
traditional RIDs requirements. Both the BD and MBD are practically equivalent 
to the WVD in the analysis and estimation of a monocomponent linear FM signal 
[5, 14], 
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6.1.5.43 Reduced interference TFDs with separable kernel filters 

A simple way to achieve a good cross-terms suppression, and at the same time pre- 
serve the signal components’ time-frequency resolution, is to design a TFD using a 
separable kernel filter expressed as 

g(r,v) = 52(f) Gi(v). (6.1.27) 

This product form allows us to formulate separate filter constraints for each variable. 

A TFD corresponding to the separable kernel filter g(r, v) is then defined, using 
Eq. (6.1.20), as 

W z (t,f) * f G 2 (/), (6.1-2S) 

where gi (£) = T ~ 1 {Gi(i/)} and G 2 (f) = ^ {g^r)}. 

IS— It T — *f 

Separable kernel filters, therefore, allow for separate convolutions of the WVD 
in both time and frequency directions. This makes the design of such TFDs easier to 
implement. Note that the (t, f) domain form of the separable kernel filter g(r.i/} is 
also a separable kernel filter 7 (f, /) = <?i(t)C?2(/) [2]. 

Two special subclasses of separable kernel filters are: 

• the Doppler-independent ( G\{u ) = 1) kernel filters, which allow for smoothing 
of the WVD in the frequency direction only, and 

• the lag-independent (#2 (r) = 1) kernel filters, which allow for the WVD smooth- 
ing in the time direction only. 

The properties of the Doppler-independent and the lag-independent kernel filters are 
studied in detail in [2, Article 5.7], 

For the WVD of a multicomponent signal, the inner interference terms [3] re- 
sulting from a nonlinear FM law of the signal components alternate in the frequency 
direction, so they can be successfully suppressed by a Doppler-independen! kernel 
filter [2, Article 5.7]. The cross-terms of the WVD oscillate in the time direction [3], 
and so they can be significantly suppressed by filtering the WVD with a time window, 
which is the inverse FT of a lag-independent kernel filter [2, Article 5.7] . To success- 
fully suppress both types of these interfering terms in the WVD of a multicomponent 
signal, a full separable kernel filter can be used. 

TheBD kernel filter [14], g(r, v) = (|r|/ cosh 2 (t)) /3 , is an example of separable 
kernels. It was shown in [14] that, for small values (close to zero) of its parameter 3. 
this TFD in general performs well for different types of signals. For small values of 
0, the BD kernel can be approximated by the lag-independent kernel filter gi(t) — 
1/ cosb 2,3 (<) [2, Article 5.7], which, when normalized, defines the MBD kernel filter 
(see Eq. (6.1.26)). The MBD achieves in particular good cross-terms suppression and 
good resolution of closely spaced components of multicomponent signals when the 
components have slowly varying IF laws, as shown in [2, Article 5.7], and illustrated 
by several examples in this chapter. The definitions and properties of other classical 
and popular quadratic TFDs are provided in [2], 
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For a more conclusive assessment of the relative worth of each TFD, a quan- 
titative comparison of their perfonnance is needed. This requires the introduction 
of specific criteria that take into account usual key attributes of TFDs (such as the 
amplitudes of signal components and cross-terms, components’ instantaneous band- 
widths, sidelobes’ amplitudes, etc.) [1], as detailed in the next section. 



6.2 Performance assessment for quadratic TFDs 

6.2.1 Visual comparison of TFDs 

Stationary signals are usually analyzed and compared in either the time or the fre- 
quency domain. The autocorrelation function in time is examined by looking at the 
position and relative amplitudes of each lobe, as well as the “correlation width” of 
the mainlobe. These characteristics are used for examination and comparison. The 
PSD in frequency is examined by looking at the position and relative amplitude of 
spectral peaks. It is desired to resolve closely located spectral peaks in the PSD. 

For nonstationary signals with a time-varying spectrum, TFSP techniques are 
needed to represent the signals in the joint time-frequency domain using an appro- 
priate choice of TFD. 

Just as some spectral estimates are better than others, some TFDs outperform 
others when used to analyze certain classes of signals [3, 15-17], For example, the 
WVD is known to be optimal for monocomponent signals with the quadratic phase 
law (linear frequency modulation (LFM)), since it achieves the best energy concen- 
tration around the signal IF law [2,3] (see Figure 6.3). The spectrogram, on the other 
hand, although still regarded as one of the most popular quadratic TFDs, results in 
an undesirable smoothing of the signal energy around its IF [3]. 

This example is just a simple illustration of the fact that choosing the right TFD 
to analyze the given signal is not straightforward, even for monocomponent signals. 
The task, then, appears to be more complex when one deals with multicomponent 
signals. 

For illustration, let us consider a multicomponent bird song signal [18], repre- 
sented in the (f , /) domain using the Bom-Jordan distribution [2], the CWD, the 
MBD, the Rihaczek distribution [10], the smoothed WVD [3, 6], the spectrogram, 
the WVD, and the Zhao-Atlas-Marks distribution [19] (Figure 6.8). 

How do we determine which of the TFDs in Figure 6.8 best represents the given 
signal? To answer this question, according to the common practice in TFSP, one 
would visually compare the eight plots and choose the one that is most appealing. 
From Figure 6.8, we can see that the MBD, the smoothed WVD, and the spectro- 
gram have “cleaner” plots (less interference, better components’ concentration) than 
the other considered TFDs. However, selecting the best time-frequency distribution 
among those three TFDs, based only on the visual comparison of their plots, is diffi- 
cult. 

The need to objectively compare the plots in Figure 6.8 requires the definition of 
a quantitative performance measure for TFDs. Some theoretical measures that deal 
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(a) Bom-Jordan 



(b) Choi-Wilhams 



(c) Modified B 



(d) Rihaczek 



(e) Smoothed Wigner- 



(f) Spectrogram 



(g) Wigner-Ville 



(h) Zhao-Atlas-Marks 



Figure 6.8. TFDs of a multicomponent bird song signal. 



essentially with signal concentration have been proposed in the literature [20-23]. 
This chapter uses objective quantitative measure criteria that take into account not 
only concentration but also resolution aspects for a practical analysis in the case of 
closely spaced components. The characteristics of TFDs that influence their resolu- 
tion, such as energy concentration, mainlobes separation, and sidelobes and cross- 
terms minimization, need to be combined to define these quantitative measure crite- 
ria, as explained in the following section. 
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FIGURE 6. 9 . Instantaneous (t = to) TFD resolution performance for a monocomponent FM 
signal. 



6.2.2 Performance criteria for TFDs 

6.2.2. 1 Monocomponent signal 



For a monocomponent FM signal, performance of its TFD is usually defined in terms 
of the energy concentration that the TFD achieves about the signal IF [4]. One desires 
to minimize sidelobe amplitude Ag(t) relative to mainlobe amplitude and to 

minimize instantaneous (mainlobe) bandwidth Vi (t) about the signal IF /;(t). 

For a given time slice of a TFD of a monocomponent signal, as illustrated in 
Figure 6.9, we may then quantify the signal TFD performance by the measure p 
expressed as 



P(t) = 



Ml 

Am {t ) 



Vj(t) 

/«(<)' 



( 6 . 2 . 1 ) 



For clarity of presentation, we limit ourselves to measuring the instantaneous band- 
width at 0. 7071 of the component normalized amplitude. Note also that in Eq. (6.2. 1) 
we have normalized the bandwidth Vi(t) with the IF /,(t). However, other possible 
normalization factors are currently being studied, such as the signal sampling fre- 
quency, or the signal half-power bandwidth [11]. 

From Eq. (6.2. 1), good performance of a TFD is characterized by a small (close 
to zero) value of its measure p. The WVD of an LFM signal with infinite duration, 
for example, has p = 0 [2, Article 7.4], 



6, 2.2.2 Multicomponent signal 

For a multicomponent FM signal, performance of its TFD involves not only the 
energy concentration the TFD attains about the respective IFs of each component, 
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Figure 6. 10. PSD estimate resolution performance for two sinusoids. 



but also the resolution, as measured by the minimum frequency separation between 
the components’ mainlobes for which their amplitudes and bandwidths are just pre- 
served. Let us define these two notions. 

6 . 2 . 2 . 2.1 Energy concentration 

By extending the concept introduced in Subsection 6.2.2. 1 for a monocomponent 
signa , we say that a TFD has best energy concentration for a given multicomponent 
signal it for each signal component it yields the smallest: 

• instantaneous bandwidth relative to the component IF (VAt)/ f At)), 

• sidelobe magnitude relative to mainlobe magnitude (|/l s It) jA M (f)'|). 

6.22.2.2 Resolution 



T a PSD eMimaB <*• *■* Of two single tones. 

fi 2inu j 2, is defined cis the minimiim p p & * • , . „ * , 

inequality holds: dtllerence f 2 - j x f or which the f ollowing 



Ji + — < f 2 



h < fi, 



( 0 . 2 . 2 ) 



as illus : tratecH n ^Figure 6 JO ^ S ° f * e first and the second sinusoid, respectively, 

accoim t rirefecToVmIs S f * tW °' COm P onent si 8 n ^ however, we also need to 

TFDs should first minimize the cross-tenns relative to ^igna^ compone^ 
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(a) WVD of the first 
component 



(b) WVD of the second (c) WVD of the sum 

component 




(f) Slice of (c) 




Figure 6.11. WVD (a~c) and its slice at the middle of the signal duration interval (d-f) for 
a two-component signal, considered individually for each contributing term and combined. 



the concentration requirement for each signal component needs to be complemented 
by the cross-term suppression requirement when evaluating the resolution perfor- 
mance of TFDs. All important signal parameters that are needed for this purpose are 
shown in Figure 6.12, which represents a time slice (t = to) of a typical quadratic 
TFD. 

In Figure 6.12, (to), f h (to), A Sl (to) and A Ml Oo) denote, respectively, the 

instantaneous bandwidth, the IF, the sidelobe amplitude, and the mainlobe ampli- 
tude of the first component at time t = to. Similarly, V %2 (to), fi 2 (to), As 2 (to), and 
Am 2 (to) represent the instantaneous bandwidth, the IF, the sidelobe amplitude, and 
the mainlobe amplitude of the second component at the same time to. Ax (to) is the 
cross-terms amplitude. 

An example of a TFD with nonresolved components is shown in Figure 6.13, 
where the signal components and the cross-term have all merged in a single lobe. 

6.2.2. 3 Resolution performance measure for quadratic TFDs 

From Eq. (6,2.2) and Figure 6.12, the frequency resolution of a TFD for a pair of 
components in a multicomponent signal may be quantified by the minimum differ- 
ence f i2 (t) - fi x (t) (/i, (t) < f i2 (t)> for which a separation measure D between the 
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Figure 6.12. Instantaneous ( t — to) TFD resolution performance for a two-component FM 
signal with resolved components. 



components’ mainlobes, centered about their respective IFs and is positive. 
We define the components’ separation measure D as 



m = 



/«(*) 



1 ( V h (t) + V ia (t) \ 

2 V/ia(*)-/ii(*)/ ' 



For a better resolution performance of quadratic TFDs: 

• the separation measure D should be maximized (close to 1), and concurrently , 

• the interfering terms (cross-terms and components’ sidelobes) should be mini- 
mized. 

To maximize D, we need to maximize the energy concentration for each component 
in a pair of signal components under analysis by minimizing their instantaneous 
bandwidths about the components IFs. 

The interfering terms, on the other hand, can be minimized by minimizing both 
the sidelobe-mainlobe amplitude ratio \As{t)/AM{t)\ (so also improving the energy 
concentration) for each component, and the cross-term— components’ mainlobe am- 
plitudes ratios \Ax{t) J AM{t)\. 

Therefore, an overall measure P of the resolution performance of a TFD for a 
pair of components in a multicomponent signal can then be expressed as [1] 



P{t) = 



As(t ) A x {t ) 1 

A M (t) Atvt{t) D(ty 
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Figure 6.13. Instantaneous (f = to) TFD resolution performance for a two-component FM 
signal with unresolved components. 



where A M (f ) , As (f ) , and Ax{t) are, respectively, the average amplitude of the com- 
ponents’ mainlobes, the average amplitude of the components’ sidelobes, and the 
cross-temi amplitude, and D(t), defined by Eq. (6.2.3), is a measure of the compo- 
nents’ separation in frequency. 

From Eq, (6.2.4) we can see that a good resolution performance of a TFD for 
a given pair of components in a multicomponent signal is characterized by a small 
(close to zero) positive value of the measure P. 

6,2.2 A Assessment of the resolution performance of TFDs for a 
two-LFM-component signal 

As an illustration of the use of the performance measure P in Eq. (6.2.4), let us 
consider a multicomponent signal si(t) of duration T = 128 given by 

s x {t) = cos (2 tt(0.15 t + 0.0004 1 2 )) + cos (27r(0.2f + 0.0004 1 2 )). (6.2.5) 

The signal s\ (f ) is represented in the time-frequency domain using a selection of 
TFDs (see Figure 6.14). In this example, we compare the TFDs’ resolution perfor- 
mance at the middle of the signal duration interval. So, for each TFD we take a slice 
at t = 64 and measure the parameters A M { 64), A§(64), Ax (64) and V)(64) (the 
average of the components’ instantaneous bandwidths). Using these, we then cal- 
culate the frequency separation measure of the components D(64), defined by Eq. 
(6.2.3), and the resolution performance measure P(64), defined by Eq. (6.2.4). The 
measurements results are recorded in Table 6.1, and the slices of the TFDs at t = 64 
are shown in Figure 6. 15. 
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(a) Bom-Jordan 



(b) Choi-Williams 

(a = 2) 



(c) Modified B 
(0 = 0 . 01 ) 
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(a) Modified B (solid), 
Bom-Jordan (dashed) 



(b) Modified B (solid), 
Choi-Williams (dashed) 



(c) Modified B (solid), 
Spectrogram (dashed) 




(d) Modified B (solid), (e) Modified B (solid), 

Wigner-Ville (dashed) Zhao-Atlas-Marks 

(dashed) 



Figure 6,15. Normalized slices of TFDs plotted in Figure 6.14 taken at a half of the time 
interval (f = 64) of signal si(t) (Eq. (6.2.5)). 



6.2.2. 5 Alternative (normalized) resolution performance measure 



An alternative to the measure P, defined by Eq. (6.2.4), is to combine | As ( t ) /Am (i ) | , 
\Ax{t)/AM{t)\, and D(t) into a sum, rather than a product, and so account for their 
effects more independently. This results in the following definition for the resolution 
performance measure [24]: 



P l (t) = 



As(t) 

A M (t) 



Ax(t) 1 

A M (t) D(ty 



(6.2.6) 



We would like to have a performance measure that is close to one for TFDs that 
perform well, and close to zero for poor-performing ones. This requires P\ to be 
normalized, for which we need to normalize each of its three contributing terms. 

Since, from Eq. (6.2.3), 0 < D(t) < 1 for resolved components, and therefore 
1 < 1/D(t) < oo cannot be normalized, we have chosen to use (1 — D(t)) instead 
of 1 /D(t) in Eq. (6.2.6). Note that we want to minimize 0 < (1 - D(t)) < 1 in 
order to achieve good frequency separation of signal components. 

On the other hand, the si delobe magnitude |A,g(t)| will always be smaller than 
the mainlobe magnitude |Ajvr(t)|, while the cross-term magnitude | Ax (t ) | can be as 
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large as twice the magnitude of the signal components [3] (we assume here normal- 
ized amplitudes). Therefore, 






< 1, and 0 < 



1 Ax ( t ) 

2 Am if) | 



By combining the three normalized quantities into a sum, and normalizing this sum, 
another resolution performance measure P 2 is obtained [24]: 



m -s(« 



, 1 Ax{t ) . 
2 A M (t) 



+ (1 ~D(t)) 



where the scaling factor 1/3 normalizes the summation. 

Smaller values of P 2 mean better resolution performance of TFDs. To make this 
measure be close to 1 for well-performing TFDs and 0 for poor-performing ones 
(TFDs with large interfering terms and poorly resolved components), we have de- 
fined the normalized instantaneous resolution performance measure Pi as [2, Article 
7.4] 



P<(t) = 1 



+ i AM 

A M (t) + 2 A M (t) 



0 < Pi(t) < 1. 



6.2.3 Performance specifications for the design of high resolution 
TFDs 

In order to design TFDs that behave optimally in the difficult case of closely spaced 

multicomponent signals, we need to revisit the constraints for TFD design, discussed 

“ ' bs ^ U0n 6 ' U ' 2 ' 3 > and relate ‘hem to the performance measure P z defined by 
cq. (o.2.y). J 

oretatinn^f TFn 1 ; 4 ’ 2 ' 3 Iists five P ro P ert ies. The first two ensure a physical inter- 
pretation of TFDs (i.e. energy conservation), while the last two directly relate to the 

Pl ' The third pr0perty ensures ** ^ IF Peak requirement is 
Sd ones pr0pertles C0Uld be further reduced t0 essentially three com- 

freauenc'Tanalv'k t001 for 3 P ractical resolution time- 

frequency analysis, a TFD must verify the following properties: 

1 6.u e ™)“ d si S“» “'“8V (Properties PI and P2 of Subsection 

2. Reveal the IF l,w(s) of a signal by its peak (s) (Propeny P3 of Subsection 

3 ' “m! 1 ! ”r“ re ?' “f “ ed by »>■ <«• 2 ! ». "> order to reduce the cross- 

n0 ' se P“fo" 
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In order to define a quadratic TFD that best meets these constraints, we could use the 
separable kernel TFD design procedure of Subsection 6.1. 5.4.3, and then vary the 
windows G \{v) and p 2 (f) until the value of the performance measure Pi is maxi- 
mum for a given class of signals. 

One outcome of such a procedure is the MBD [5], whose kernel filter is defined 
by Eq. (6.1.26). The kernel filter of the MBD was chosen in the ambiguity domain to 
be a 2D low-pass function centered around the origin, with sharp cutoff edges. In this 
way, the kernel retains as much of the auto-terms as possible, while attenuating cross- 
terms. The amounts of auto-terms and cross-terms kept and filtered out are functions 
of the volume underneath the kernel filter g{r,u). This volume can be changed by 
varying the parameter /3. 

The MBD and its original version the BD, meet the properties listed above, al- 
though, like the spectrogram, they do not satisfy the marginals [5, 14,25]. 



6.3 Selecting the best TFD for a given practical 
application 

6. 3. 1 Assessment of performance for selection of the best TFD 

6.3. 1.1 Assessment and selection procedure 

The procedure for selecting the optimal time-frequency representation for a given 

multicomponent signal consists of the following steps: 

1 . Define a set of criteria for comparison of TFDs: The criteria must be related 
to the information we seek from a TFD (e.g., the number of signal components, 
their relative amplitudes, components’ frequency modulation laws, etc.). A set of 
such criteria is defined in Subsection 6.2.2 for both mono- and multicomponent 
signals. 

2. Define a quantitative measure for evaluating TFDs performance based on these 
criteria: The measure Pi (Eq. (6.2.9)), for example, has been defined for eval- 
uating the resolution performance of quadratic TFDs based on the comparison 
criteria for FM signals described in Subsection 6.2.2. 

3. Optimize TFDs to match the comparison criteria as closely as possible: The 
optimization procedure for a TFD with parameter a (e.g., MBD with parameter 
/?, or CWD with parameter <r) can be done as follows. First, we choose the initial 
value for the parameter a (a value closest to its lower bound) and calculate the 
TFD of a given signal. For each time instant in the time interval of interest (time 
instants over which we want to compare the performance of different TFDs), we 
take a slice of the TFD and measure its instantaneous performance (e.g., use Pi 
defined by Eq. (6.2.9)). The average of all instantaneous measures defines the 
interval performance measure of the TFD for the given value of a. 

The procedure is repeated for the next value of the parameter a. The increment in 
a should be neither too small (long computation time) nor too large (too “coarse” 
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f 



optimization results). The optima] value of the parameter a is that value of a 
that maximizes the TFD interval performance measure. This maximum value of 
the TFD interval performance measure is called the TFD optimal performance 
measure. 

4. Select the best TFD for the given signal'. When all TFDs used to represent the 
given signal in the joint (f, /) domain are optimized, the TFD with the largest 
value of its optimal performance measure is selected as the best TFD for the 
analyzed signal. 

6.3.1 .2 Resolution performance comparison of TFDs for a two-LFM-component 
signal in additive white Gaussian noise 

The following example illustrates how to use the above procedure for selecting the 
optimal TFD for the two-component signal si(t) (Eq. (6.2.5)) embedded in noise, 
i.e., 

s 2 (t) = cos (2rr(0.15f + 0.0004 1 2 )) + cos (2tt(0.2 t + 0.0004 1 2 )) + »(<), 

(6.3.1) 

where n(t) is additive white Gaussian noise, with a signal -to-noise ratio (SNR) of 
10 dB. 

The signal s 2 (t ), of duration T = 128, is analysed in the (t, /) domain using the 
Bom-Jordan distribution, the CWD, the MBD, the spectrogram, the WVD, and the 
Zhao-Atlas-Marks distribution. 

To find the TFD that best resolves the two LFM components of the signal s 2 (t), 
we first optimize each of the considered TFDs, as defined in step 3 of the procedure 
in Subsection 6.3. 1.1 (note that the WVD and the Born-Jordan distribution (a [2] 
equals 1/2, by convention) have no kernel filter parameters, and hence need no op- 
timizing). The resolution performance of TFDs is compared based on the criteria 
defined in Subsection 6.2.2, using the comparison measure Pi (Eq. (6.2.9)). The op- 
timal resolution performance measure, P opt , is found for each of the TFDs. Table 6.2 
contains the results of the optimization process. It shows that the optimal TFD for 
signal s 2 (t) is the MBD with parameter f3 = 0.002, as it has the largest value of 
P op t. The time-frequency plots of the optimized TFDs ar e shown in Figure 6.16. 



Born-Jor dan 
Choi-Williams 
Modified B 



Optimal kernel filter parameter 



a = 0.5 



13 = 0.002 

Specteogram Bartlett window, length 47 0.8448 

Wigner-Ville N/A 0.6694 

Zhao-Atlas-Marks a = 8 0.6550 

Table 6.2. Optimization results for TFDs of signal S 2 {t) defined by Eq. (6.3.1). 



0.7542 

0.7976 

0.8605 

0.8448 
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(a) Born-Jordan 



(b) Choi-Williams 
(cr = 0.5) 



(c) Modified B 
(J3 = 0.002) 




Figure 6. 16. Optimised TFDs of signal S 2 (t) defined by Eq. (6.3.1). 



As explained in Subsections 6.1. 4.2.3 and 6.2.3, for a practical time-frequency 
analysis, an important property that a TFD should satisfy is to accurately reveal the IF 
laws of signal components by its peaks. In Figure 6.17 we compare the true IF laws 
of the two signal components with those measured from the peaks of the optimized 
MBD (the best-performing TFD for s 2 (t)), and the optimized spectrogram (second 
best TFD for S 2 (t)). 

The quality of the components’ IF estimates is measured using the mean-squared 
error (MSE). As indicated by the MSE values, recorded in Figure 6.17 next to each 
of the two components, the MBD provides more accurate estimates (smaller MSEs) 
than the spectrogram for both components of the noisy two-component signal s 2 (t). 

6.3.1 .3 Resolution performance comparison of TFDs for a signal with a linear and 

a nonlinear FM component embedded in additive white Gaussian noise 

The optimal TFD methodology is now illustrated for a two-component signal s 3 (t) 
comprising a quadratic FM component whose frequency varies from 0. 1 Hz to 0.3 Hz 
(the sampling frequency is f a = 1 Hz) over the time interval t £ [1,128], and 
a component that linearly decreases in frequency from 0.4 Hz to 0.275 Hz over the 
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(a) Modified B: MSEi = 2,09 x (b) Spectrogram: MSEi = 3.8(1 >- 
1CT 6 , MSEi = 2.89 x 10 ~ 6 11T 6 , MSDa = 4.79 xl0~ 6 



Figure 6.17. Comparison of the measured (dashed) and true (solid) IF laws of the two linear 
FM components of signal Si(t), defined by Eq. (6.3.1), for the optimized MBD (a) and the 
optimized spectrogram (b) of 82 (f). 



same time interval but has nonzero rectangular amplitude values only for t G [S3. 96' ■ 
A diagram illustrating the IF laws of the signal s 3 (t) is shown in Figure 6.18. 

We measure the optimized performance P op t of the Bom-Jordan, Choi-Williams. 
Modified B, spectrogram, Wigner-Ville, and Zhao-Atlas-Marks distributions, when 
these TFDs represent s 3 (f) in the joint time-frequency domain. Table 6.3 lists the 
TFDs’ optimal performance measures and the corresponding optimal values of their 
kernel filter parameters. It shows that the signal best-performing TFD is the spectro- 
gram (P op t = 0.84), closely followed by the MBD (P op t = 0.8361). 

If we repeat this experiment, but embed the signal s 3 (t) in additive white Gaus- 
sian noise for two different values of SNR, then: 




Figure 6 . 1 8 . Diagram illustrating the IF laws of signal S 3 (f ) . 
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TFD 


Optimal kernel filter parameter 


Popt 


Bom-Jordan 


N/A 


0.7114 


Choi- Williams 


a = 0.7 


0.7660 


Modified B 


/3 = 2.6 x 10" 3 


0.8361 


Spectrogram 


Hamming window, length 53 


0.8400 


Wigner-Ville 


N/A 


0.6232 


Zhao-Atlas-Marks 


a — 7 


0.6359 



Table 6.3. Optimization results for TFDs of signal S 3 (f ). 



TFD 


Optimal kernel filter parameter 


Popt 


Bom-Jordan 


N/A 


0.7073 


Choi-Williams 


a = 0.5 


0.7414 


Modified B 


f3 = 7 x 10" 3 


0.8316 


Spectrogram 


Bartlett window, length 57 


0.8401 


Wigner-Ville 


N/A 


0.5791 


Zhao-Atlas-Marks 


a = 7 


0.6471 



Table 6.4. Optimization results for TFDs of signal S3 (t) in 5 dB additive white Gaussian 
noise. 



TFD 


Optimal kernel filter parameter 


Popt 


Born-Jordan 


N/A 


0.6548 


Choi-Williams 


<7 = 45 


0.5811 


Modified B 


(3 = 2.6 x 10~ 3 


0.8219 


Spectrogram 


Hanning window, length 47 


0.7955 


Wigner-Ville 


N/A 


0.5647 


Zhao-Atlas-Marks 


a = 2 


0.5703 



Table 6.5. Optimization results for TFDs of signal S 3 (t) in 0 dB additive white Gaussian 
noise. 



1. For SNR = 0 dB, the optimal TFD is found to be the MBD, whose optimal res- 
olution performance measure is P opt = 0.8219. The spectrogram is the second 
best TFD with P opt = 0.7955. Figure 6. 19 shows the plots of the TFDs of signal 
53 (f) in 0 dB additive white Gaussian noise. 

2. For SNR = 5 dB, the optimal TFD is the spectrogram ( P opt = 0.8401), while 
the MBD is the second best TFD ( P op t = 0.8316). 

The performance of other TFDs for 53 (f) in 5 and 0 dB noise is shown in Tables 6.4 
and 6.5, respectively. 
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(a) Born-Jordan (b) Choi-Williams (c) Modified B 

{a = 45) (/3 = 0.0026) 




(d) Spectrogram (e) Wigner-Ville (f) Zhao-Atlas-Mnrks 

(Hann, L = 47) = 21 



Figu r e 6 19, Optimized TFDs of signal s 3 (t), with the components’ IF laws defined in 
Figure 6.18, embedded in 0 dB additive white Gaussian noise. 



me optimized spectrogram and the optimized MBD of signal s 3 (t) perform (oul 
and m noise) similarly, with both TFDs clearly outperforming all other consider 

i u S ; Ua : i0nSwhen two or more TFDs have similar performance for a given 
signal, the final choice of the signal’s optimal TFD should be based on application- 
s'^ ?r n ‘ S ' S ° * at th£ TFD * at more closel y ^ets these constraints is 
ht MRO fl° Ver ° th£r “° ptimal ” TFD <»' For sample, we have seen that 

I ff) How? f Sp£Ct ^ gram haVe similar resolution performance for the signal 
d Ld in Th’ alT t0 ^ res ° Iuti0n requirement, the peak IF criterion 

Lidhertrif ^ r U } 1S Considered ’ the MBD outperforms the spectrogram, 
TFl ! ! X J 35 Ae ° Ptimal TFD f0r (*)• The issue of selecting the optimal 

example, to dlSCUSSed ’ and illustrate d on a real-life signal 



6 ' 3 ' U with Lt r rf Utl ° n Perf ° manCe ° fTFDs f0r a ^-component signal 
h different frequency separation between the components 

We now compare the resolution performance of TFn, , • , . u , , 

quency components, with different sep^rions ill ^ Wlth ^ fK ' 

parations between its consecutive components 
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in the time-frequency plane. The measure Pi, defined by Eq. (6.2.9), evaluates the 
TFD performance for a given pair of the signal components. 

For this test, it suffices to consider a signal S 4 (t), which is a sum of three sinu- 
soids: fi = 0.1 Hz, f'l = 0.2 Hz, and / 3 = 0.4 Hz, all three defined overt e [1,128]. 
We define P M to be the resolution performance measure for the first pair of compo- 
nents {/i, / 2 } (closer components), and P i2 the resolution performance measure of 
the second pair {/ 2 , fs,} (more separated components). Table 6.6 shows the values of 
P il (64) and P i2 (64) for the Born-Jordan, Choi- Williams, Modified B, spectrogram, 
and Zhao-Atlas-Marks distributions. We do not consider the WVD in this example 
since the cross-terms significantly degrade its performance whenever there are more 
than two components in the signal. 

Table 6.6 indicates that the MBD is the best-performing TFD for the pair of closer 
components, followed by the spectrogram, Choi-Williams, Bom-Jordan, and Zhao- 
Atlas-Marks distributions. The same ranking of the considered TFDs is obtained for 
the pair of more separated components. 

Thus, for signals with more than two components, we only need to determine a 
TFD’s resolution performance measure for the pair of closest components. This value 
of the performance measure provides an indication of the TFD’s resolution perfor- 
mance, with respect to other TFDs, for all other pairs ol consecutive (in frequency) 
signal components. 

6 . 3 . 1 . 5 Testing the resolution performance of TFDs for a signal with components 

of different amplitudes 

We have shown that the spectrogram and the MBD are, based on the resolution per- 
formance measure P u two best-performing TFDs for representing the signal s^f) 
consisting of three sinusoids. We now assess how the performance of these TFDs 
varies when the components 5 amplitudes of a two-sinusoid signal 55 (t) ■ (a) are equa , 
and (b) the second component amplitude is half of the first component amplitude (in 
time); both cases are studied for the following three separations between the two 
sinusoids: A / = 0.2 Hz, A / = 0.1 Hz, and A/ — O.Oo Hz. 



TFD 


(64) 


Pi, (64) 


Born-Jordan 


0.8887 


0.9029 


Choi-Williams 


0.9387 


0.9580 


Modified B 


0.9566 


0,9675 


Spectrogram 


0.9520 


0.9631 


Zhao-Atlas-Marks 


0.8054 


0.8128 



(P i2 (64)) signal components. 
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A 2 /A 1 








L =31 


Rectangular 


0.2 


1 










1/2 


0.8235 








0.1 


1 


0.8952 


0.8643 
















0.05 




0.8279 






0.8010 








SpifStp 


0.5116 


Hamming 


0.2 


1 


0.9682 


0.9646 


0.9570 


0.9316 


1/2 


0.9644 


0.9589 


0.9519 


0.9282 


0.1 


1 


0.9576 


0.9545 


0.9381 


0.9151 


1/2 


0.9534 


0.9495 


0.9334 


0.9129 


0.05 


1 


0.9403 


0.9163 


0.8682 


0 


1/2 


0.9351 


0.9175 


0.8782 


0.8704 



Table 6.7. Assessment of the resolution performance (at t = 64) for the spectrogram of a 
two-sinusoid signal ss(t) for different ratios of the components’ mainlobe amplitudes A'i/At 
and different components’ frequency separations A/. The spectrogram is calculated using the 
rectangular and Hamming windows of length L. 



Table 6.7 summarizes the performance results (at the time instant t = 64) ob- 
tained for the signal spectrogram with the rectangular and Hamming windows (of 
four different lengths). We observe that the resolution performance of the spectro- 
gram with the Hamming window continuously decreases as the window length de- 
creases from L = 127 to L = 31 for both equal and nonequai amplitudes of the 
components for all three frequency separations between these components. The dif- 
ference in the P, values becomes larger as the components get closer in frequency, 
and so their resolution becomes more critical. The spectrogram with the rectangu- 
lar window does not show this trend of the continuous decrease in its resolution 
performance, since, unlike the Hamming window case, its optimal window length 
{L op t = 107) is smaller than the signal duration T = 128. 

Table 6.8 gives the performance results (at the time instant t = 64) for the MOD 
of signal s 5 (t) for a range of its kernel filter parameter values. It indicates that for 
both the equal and nonequal amplitude components the resolution performance of 
the MBD continuously decreases as the parameter /? increases, regardless of the 
components’ frequency separation, This is due to the fact that, for large /?, cross- 
terms are less attenuated by the more-spread kernel filter in the ambiguity domain, 
and the large values of the cross-term amplitude Ax in the (£, /) domain reduce the 
resolution performance measure P*. 

This test illustrates that the measure Pi provides a good assessment of TFDs’ res- 
olution performance for signals with different components’ amplitudes and different 
time- varying components’ frequency separations. 
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1 

II 

<1 


AijA\ 


p = 10~ 4 


P = io ~ 3 


p = 10 -2 


p = w~ 1 


0.2 


1 


0.9700 


0.9211 


0.8757 


0.8683 


1/2 


0.9664 


0.9054 


0.8365 


0.8263 


0.1 


1 


0.9611 


0.9236 


0.8933 


0.8813 


1/2 


0.9568 


0.9065 


0.8622 


0.8590 


0.05 


1 


0.9436 


0.8938 


0.8520 


0.7011 


1/2 


0.9412 


0.8822 


0.8199 


0.4610 



Table 6.8. Assessment of the resolution performance (at t = 64) for the MBD (for four 
values of its kernel filter parameter /3) of a two-sinusoid signal ss(t) for different ratios of the 
components’ mainlobe amplitudes A^/Ai and different components’ frequency separations 

A/. 



6.3.2 Selecting the optimal TFDfor real-life signals under given 
constraints 

As mentioned earlier, different TFDs display information in the time-frequency plane 
with a varying amount of detail and “accuracy.” In an application, selecting a TFD 
that does this in the optimal way is a critical factor when applying time-frequency 
analysis to real-life signals. 

Based on the concept of multiple view TFDs [26], we define in this subsection 
the methodology for the selection of the optimal TFD, under application-specific 
constraints, for real-life signals in time-frequency regions where signal features of 
interest are located. 

6.3.2. 1 Methodology 

The methodology consists of the following steps. 

1 . Represent the signal in the (t, /) domain using the WVD, the spectrogram, and 
the MBD. These three TFDs will provide us with indications of the main signal 
features in the time-frequency plane: the number of components, their relative 
amplitudes, the components’ durations and bandwidths (features obtained from 
the MBD and the spectrogram), as well as the cross-terms locations (obtained 
from the signal WVD). 

The WVD is essentially a parameter-free TFD, while the key parameters of the 
spectrogram and MBD (the window type and length, and the kernel filter pa- 
rameter P, respectively) need to be initialized. For the spectrogram, we use the 
Hanning window as the initial window type, and from several window lengths 
(short, medium, and long — whose values depend on the signal time duration) 
select as the initial window length the one that results in the “cleanest” time- 
frequency plot. This window length is also used as the initial effective window 
length [3] for all other TFDs considered in the application. In the case of MBD, 
the initial value for P is set to 0.001. The choice of these initial parameter values 
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is based on the fact that the spectrogram with the Hanning window and MBD 
with (3 < 0.001 have been found to perform well for the majority of signals 
whose TFDs’ resolution performance we have analyzed in [1,2,27,28]. 

2. Using these three TFDs, identify regions in the (1, /) domain where application- 
specific signal features (e.g., closely-spaced components, crossing components, 
sudden changes in component IF laws) are located. We call such regions the 
regions of interest (ROIs), The ROIs are rectangles in the time-frequency plane 
whose dimensions are the “time of interest” and the “frequency of interest.” The 
time (resp. frequency) of interest is the time interval (resp. frequency band) of 
the TFD where the specific signal features are located in time (resp. frequency). 

3. For each selected ROI, optimize different TFDs under application-given con- 
straints. 

4. Select as the signal optimal TFD, the TFD which among other considered opti- 
mised TFDs best satisfies each of the given constraints in all observed ROIs. 

6.3. 2.2 Optimal TFD selection for the unbiased, efficient, multicomponent IF 
estimation of the Noisy Minor song signal 

In this section we illustrate the use of our methodology for the “Noisy Minor” ( Mono - 
rina melanocephala ) Australian bird song signal [18], The constraint that needs to 
be satisfied is the unbiased, efficient, multicomponent IF estimation in (t, f) regions 
where the signal components are closely-spaced. 

Figure 6.20(a) shows the signal time-domain plot, which indicates how its ampli- 
tude varies with time, and Figure 6.20(b) (the signal PSD) shows which frequencies 
(with which magnitudes) are present in the signal. However, neither of the two can 
provide us with information on the signal internal structure (the number of signal 




(a) Time -domain 



(b) Frequency-domain 



Figure 6.20. Time-domain and frequency-domain representations of the Noisy Minor song 
signal. 
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components, their relative amplitudes, IF laws, time intervals and frequency bands 
these components occupy, etc.). 

To have a more complete “picture” of this signal, we analyze it in the [t , /) 
domain. This requires identifying a TFD that is optimal for the signal. By employing 
the steps of the above-defined methodology, we find in this example the regionally 
optimal TFD of the Noisy Minor song signal for efficient and unbiased estimations 
of the signal components’ IFs from the TFD’s peaks [2, 3, 6] in the (i, /) regions 
where these components form closely spaced pairs. 

We start by representing the signal in the (t , /) plane using the WVD, MBD (/? = 
0.001), and the spectrogram with Hanning window of length 511. In selecting the 
optimal window length for the spectrogram, other lengths (127, 255, 1023) were also 
tested, but length 511 resulted in the most clear (t , /) plot of the signal spectrogram. 
The effective window length for the WVD and the MBD is also initially set to 511. 
By comparing the plots of these TFDs, several components (dominant ridges) are 
identified in the signal, and four ROI are defined where these components are closely 
spaced to each other in the time-frequency plane (see Figure 6.21). 

In order to meet the accurate IFs estimation constraint, we first need to optimize 
the resolution performance of TFDs in the selected ROIs by applying the proce- 
dure defined in Subsection 6.3. 1.1. Table 6.9 lists the TFDs’ resolution optimiza- 
tion results. It shows that the MBD with (3 = 5 x 1 0 ” 5 , having the largest optimal 
performance measure P opt (the regional performance measure corresponding to the 
optimal value of the TFD kernel filter parameter) among the eight considered TFDs, 
achieves the best resolution of the signal components in ROIi . It slightly outper- 
forms the smoothed WVD (with the Hamming window of length 415 chosen as 




Figure 6.21. ROIs for the Noisy Minor song signal’s time-frequency distributions. 
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the smoothing window) and the spectrogram (with the Hanning window of length 
447). Other TFDs, even after being optimized in ROIi, still do not achieve as good 
resolution performance as the MBD does. The Rihaczek distribution and WVD in 
particular perform poorly, as indicated by their P opt values of 0.5825 and 0.4984, 
respectively. 

Figure 6.8 shows the ROIj of the optimized TFDs of the Noisy Minor song sig- 
nal. Note that P opt ranks the TFDs in accordance with our visual impressions of 
their (f, /) plots; i.e., the clearer the plot is (good components’ concentration and 
less interference present in the TFD), the larger P opt is obtained. 

Similar analyses can be done for ROI2, ROI3 and ROI4. What is interesting 
to observe is that the spectrogram has the best resolution performance in each of 
ROI2, ROI3, and ROI4. The spectrogram’s performance closely approaches that of 
the MBD in ROIi too. The three best-performing TFDs for all ROIs considered in 
this example, with very similar resolution performances, are the MBD, the smoothed 
WVD, and the spectrogram. The fact that more than one TFD perform well for a 
certain signal can be of benefit to the signal analyst, giving him/her more freedom in 
selecting the signal optimal TFD under application-specific constraints. 

Having dealt with the requirement for a good (t, f) resolution of closely spaced 
components of the Noisy Minor song signal, we now select the signal TFD that al- 
lows for the most accurate components’ IFs estimation (from the peaks of the TFD’s 
dominant ridges) in each of the observed ROIs. Among the three best-performing 
TFDs, only the MBD provides unbiased and efficient multicomponent FM signals’ 
IF estimates [5], Therefore, under the given accurate multicomponent IF estimation 
constraint, we select the Modified B distribution as the optimal TFD of the Noisy 
Minor song signal, for all four ROIs, The optimal value of the MBD kernel filter 
parameter (3 slightly varies across the ROIs, so that in each of these regions the best 
resolution of the signal components is achieved; as required by the robust multicom- 
ponent IF estimation techniques [2]. 



6.4 Discussion and conclusions 



This chapter brings together, in a heuristic way, three key recent developments that 
are fundamental to a better understanding and use of time-frequency signal analysis 
tools. 

The first development consists of a concise, updated presentation of the core 
concepts of TFSP, including a simple way of designing high resolution TFDs for 
real-life practical applications. 

The second development consists of using a quantitative measure of goodness 
for TFDs to compare the resolution performance of time-frequency distributions for 
multicomponent signals analysis. This result fills an obvious need for the practitioner 
in that, until recently, the comparison of the resolution performance of TFDs was 
primarily based on a visual impression of the TFDs’ plots. The introduction of the 
resolution performance measure for quadratic TFDs has led to an improvement in 
the design of tools for high resolution time-frequency analysis of multicomponent 
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1 
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1 


0.7922 
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Table 6.9. The optimal resolution performance measures P op t and the corresponding optimal kernel filter parameter values of 
TFDs of the Noisy Minor song signal for the four time-frequency ROIs defined in Figure 6*21. 
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signals, by removing unnecessary limitations in the way desirable properties of TFDs 
were previously chosen. The MBD, which outperforms existing TFDs in terms of 
(t, /) resolution for signals with closely spaced components, was introduced in this 
way. 

Thirdly, a methodology for selecting an optimal TFD for a given real-life sig- 
nal under application-specific constraints has been defined. The use of the method- 
ology in practice was illustrated on an Australian bird song signal, for which the 
optimal TFD for an accurate multicomponent IF estimation was found in several 
time-frequency ROI. 

The results presented in this chapter are important for the field of TFSP. They 
not only allow for the optimal selection of a specific application-dependent TFD, but 
also open the way for further research in developing high resolution digital signal 
processing tools for nonstationary signals, by providing a measure of quality of TFDs 
and removing often unnecessary TFD design limitations. 
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